Low-temperature proximity effect in clean metals 
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Theories of proximity effect in layered superconductor-normal metal (SN) structures usually deal 
with a hypothetic normal metal with no interaction between electrons and with finite temperatures 
often close to the superconductor critical temperature. We present an asymptotic solution of the 
Gor'kov equations in the opposite low-temperature limit for a clean normal metal with a repulsive 
interaction between electrons. The order parameter in the metal exhibits a power-law decay, A(x) oc 
£/x, as a function of the distance from the SN boundary, x, with a proximity length £ strongly 
depending on the repulsive interaction. 

PACS numbers: 71.38.-k, 74.40.+k, 72.15.Jf, 74.72.-h, 74.25.Fy 



In recent years investigations of different SN struc- 
tures [l[ have gone through a vigorous revival. In par- 
ticular, superconductor /ferromagnet structures 0, Q, 
cuprate SNS jun ctions [l| @| , and mesoscopic SN struc- 
tures @, H, H,H3| have been experimentally studied and 
addressed theoretically. 

The superconducting order parameter penetrates into 
a bulk normal metal across the SN boundary. The micro- 
scopic theory of this proximity effect has been developed 
at finite temperatures (for reviews see [ill, EH), using 
the Eilenberger formalism [l3| and the semiclassical Us- 
adel approximation [3] for solving the Gor'kov equations 
15 1 close to and below T c , in particular in the dirty limit 
16| . The clean case has been studied by Falk [17| for 
a hypothetic normal metal with no interaction between 
electrons, and by Maksimov and Potapenko [l|| for a 
"normal" metal with a weak attractive interaction and 
reduced critical temperature. An exponential decay of 
A(x) oc exp(— x/£), as a function of the distance from 
the SN boundary, has been found with the proximity 
length £ cx 1 /T in the clean limit [l8[ , and (oc 1/ y/T in 
the dirty limit [l], HH, 03 • To the best of our knowledge a 
role of the repulsive interaction in the normal metal has 
not been addressed in the clean case, while its effect on 
the proximity length and the Josephson current has been 
found marginal in the dirty case [1 €31 ] - 

Here we present an asymptotic solution of the Gor'kov 
equations for the SN boundary between a bulk supercon- 
ductor and a bulk clean normal metal with the repulsion 
between electrons at low temperatures, Fig.l. 

It is convenient to fourier-transform the Mat- 
subara normal, Q u (r — r' ,x,x'), and anomalous, 
•F<j"(r — r',x, x'), Green's functions (GFs) along the 
boundary, r = {y, z} with the wave vector k = {k y , k z }, 



(2tt) 2 J dkQ Uik (x,x')exp[ik- (r 
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dkJ-* k (x, x') exp[ik ■ (r — r')]. 



tion for the Matsubara operators as 
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— I « 2 + ) G u , k (x,x')+A(x)F+ k (x,x') = S(x-x'), 
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~L ( a * 2 + &) - A*(x)g u , k (x,x') = o, 



(2) 
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Here m, /cf are the electron effective mass and the Fermi- 
momentum, respectively, which are taken the same in the 
superconductor and in the normal metal for mathemat- 
ical transparency, a 2 — 2m(iuj — £fe) with £fc = [k — 
k 2 F )/2m and a = sign(u) yjm(p + + iy / m(p — 
p = +y/uP + £f, uj = 2nT(n + 1/2) is the Matsubara 
frequency (n = 0, ±1,±2, ...). We use h = ks = 1 here 
and below, and !F_* k (x, x') — T U) k{x,x'). The super- 
conducting order parameter 

A{x) = ~V(x)T(2tt)- 2 [ f&fuAx, x), (5) 
is a solution of an integral equation, 

/oo 
dx'g m uk (x,x')A(x')g M , k (x,x'), 
-oo 

where V(x) is the electron-electron contact interaction, 
which is negative (attractive) in the superconductor at 
x < 0, and positive (repulsive) in the normal metal, 
V(x) = V c > 0, at x > 0, Fig.l. The GF of a bulk 
normal metal is 



The Gor'kov equations are derived using equations of mo- 



Gu?k( x > x ') = — exp(ia\x - x'\). 



(7) 
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FIG. 1: The order-parameter profile near the boundary be- 
tween a superconductor (S) and a normal metal (N) with the 
repulsive electron-electron interaction. 




FIG. 2: Four domains with different solutions for GFs with 
the step-like order parameter A(x) = A s 0(— x). 



The Gor'kov equations are supplemented by the bound- 
ary conditions, A(— oo) = A s , A(oo) = 0, and all GFs 
should be continuous with respect to x, x' together with 
their first derivatives at the boundary x = or x' = 0. 

To solve the integro-differential system of equations 
[T]|n] let us assume that the repulsive interaction in the 
normal metal significantly reduces the order parameter, 
so that the latter is small, A(x) <C A s , far away from 
the boundary at x 3> £, where the proximity length 
£ is small compared with the superconductor coherence 
length, £ s , £ <C £ s - That allows us to use in Eq.© a so- 
lution, Guj,k(x',x) w Guk{x',x) of the Gor'kov equations 
[T]|4]with a step-like order parameter A(x) = A s 0(— x), 
where 0(x) is the Heaviside step function. This solu- 
tion can be readily obtained by matching GFs and their 
derivatives at the boundaries between 4 domains: a su- 
perconductor domain S, where both arguments are neg- 
ative x, x' < 0, two mixed domains M with x < < x' 
or x 1 < < x, and a normal domain N , where x, x' > 0, 
Fig.2. 

In the normal domain the solution is found as (see also 

[am) 



g^ k {x,x')+Ae ia ^ +x '\ 



T+^ N (x.x') = Be i{ - ax '- a * x \ 
and in one of the mixed domains, x' < < 

2mA s t(ax+b * x ') 
+ (6* 2 -a 2 ) 



(8) 



(9) 



(10) 



F+£ M (x,x') = c*e t{h * x '- a ' x) + D*e' l{bx ' +a ' x \ (11) 

where A,B,C and D are constants and b 2 = 2m(ie — 
36 > 0, e = \J ijS 1 + A 2 . The constants are found from 



SZiUx, o) = gffi (x, o), (x, o) = (*> o) ^d 

iromJLg^(x,x>) = '(x,x% ^T^ N {x,x') = 

£r^t tM (x,x') at x' = 0. 



In particular we find 

Tji ( 2a((a* + b*)(e + us) + (a* 



A = 



b){e-u)) 



\a + b\ 2 {e + cu) + \a-b*\ 2 (e-Lu) 



and 



B = 2mA, 



b + b* 



\a + b\ 2 (e + uj) + \a~b*\ 2 (e-Lu)' 



(12) 



(13) 



Integrating over the normal region, x 1 > in Eq.© 
one can keep only the first "normal" term of Eq.© at 
sufficiently large x > £, while integrating over the super- 
conductor region, x' < 0, one can use directly Eq.© and 
the definition of A (a;), Eq.© to obtain 



dx'A(x')e- 2 ^- x '^ a 



-2xQa 



(14) 



where A = A(x)/A s is the reduced order parameter in 
the normal metal, and u = B/mA s . 

At finite temperatures and x 3> vf/(2ttT) the main 
contribution to the Matsubara sum in Eq. (fl"4|) comes 
from the n — term, so that the order parameter 
has the conventional exponential asymptotic, A(x) oc 
exp(— 2ttTx/vf), where Vf — hp/m is the Fermi veloc- 
ity. At sufficiently low and zero temperatures one has 
x -C vf/(2itT) for any size of the normal region, so that 
the exponential asymptotic is replaced by some power 
decay [19]. To find the power we replace the Matsubara 
sum by an integral over u> and integrating over momen- 
tum and frequency in Eq. (114|) obtain an integral equation 
for the reduced order parameter, A(x) = A(x)/A s , 



dx 



, A(x') 



A(x) 



(15) 
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FIG. 3: SN proximity length £ (in units of the superconductor 
coherence length) as a function of the repulsive pseudopoten- 
tial (U c (solid line) and of the attractive potential (dashed 
line) . 



where /i c = V c mkF/^ 2 is the repulsion pseudopoten- 
tial, and £ s — i^f/2A s . One can satisfy Eq. (fT5|) with a 
solution decaying as inverse distance from the boundary, 



A(x) = -i/x. 

Substituting Eq.([T6]) into Eq.Jl5]) yields 

Mr 



i + fJLc j™dt/t\t-i\- 



(16) 



(17) 



where £ = £/£ s is the dimensionless proximity length. 

The logarithmic divergency of the integral in Eq. (|17|) 
is an artifact of the step-function approximation, Eqs.© 
[9]). We cut the divergency by excluding small regions, 
t < t m i n <C 1 and \t — 1| < t m i n from the integral, where 
the step- function approximation fails because A (a;) be- 
comes relatively large. We chose t m i n — £, since £ is the 
only dimensionless parameter in the normal region which 
corresponds to the cutoff of the integral over x' in Eq. (jX5jl 
at a cutoff length x' > l m in = proportional to x. The 
result is a transcendental equation for £, 



l-3/x c ln(0 



(18) 



The proximity length is shown in Fig. 3 as a function of 
the repulsion fi c . It strongly depends on the repulsion 
with a maximum value £ ma z ~ 0.22 at /i c = oo [20j. The 
magnitude of £ is small at any /i c , which justifies our 
step- function approximation, Eqs.([51 O for solving the 
problem. 

To verify the self-consistency of the approximation one 
can estimate the correction, SGui,k(x,x'), to GF, Eq.© 
due to the finite order-parameter in the normal region, 



at large x, x' > 0, 

/>oo 

sg u ,kM = / dx'gW( x , x ")A(x^+t N (x", x '), 

Jo 

(19) 

Using Eq.© and Eq. p^j) one obtains 8Q u ^{ X ,x') cx 
£ln(x/l m i n ), which is small as £ln(l/£) < 1 at any x 
with our choice of the cutoff length, l m in — t x - An- 
other possible choice of the cutoff length l m i n = £, 
does not change the order parameter in a wide region 
£ s < x < ^ s /^ because the singularity is logarithmically 
weak. We note that the order parameter becomes so 
small at very large distance from the boundary, that the 
BCS mean-field approximation used in Eqs. (jTJ- [4j) may 
break down plj . 

While the order parameter saturates at a small value 
for any fixed x in the normal region, the pair wave func- 
tion, !F^ k (x,x'), decreases as l//i c at large repulsion. At 
first glance the existence of finite pair correlations in the 
repulsive normal metal looks surprising because it leads 
to a finite increase of the repulsive potential energy of the 
whole SN system. In particular, an individual real-space 
pair of electrons bound by some attractive potential on 
one side of the boundary would stay at an infinite dis- 
tance from the boundary with the repulsive interaction 
(zero proximity effect). However, the electron-density ho- 
mogeneity in metals creates a quantum pressure on pairs 
pushing them across the boundary. The reason for the 
failure of the pair wave function to die off is that there 
is simply no pair-breaking mechanism for disrupting any 
correlation that drift across the boundary as in the case 
of a hypothetical normal metal with no interaction [l7l ]. 
With the power law order parameter, A(x) oc 1/x, an 
increase of the potential energy, proportional to A 2 (x) 
is finite, when it is integrated over the whole normal re- 
gion. Compared with the step-like order parameter this 
increase is compensated by a lowering of the kinetic en- 
ergy near the boundary. 

It is instrumental to compare the SN proximity effect in 
the repulsive normal metal with the effect in a " normal" 
metal with a small attractive potential between electrons 
at low but finite temperatures above the transition tem- 
perature, T cn (<C T c ), of the "normal" metal. Solving 
Eq. (fT5)) with a negative \i c yields a positive order param- 
eter 



A ott (a;) = iattjx 



with the proximity length, 



£,att 
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l + 3|Mc|ln(Cat 



(20) 



(21) 



shown in Fig. 3 as the dashed line. With increasing attrac- 
tion the proximity length increases linearly as in the case 
of the repulsion, but much faster saturating at ^ aU ~ 1.2 
for the values of the interaction which are well beyond 
the step-function (and BCS) approximations. The power 
law decay of the order parameter, Eq. (j2"D|) , holds for the 



4 



region x <C vf/(27tT), but different from the repulsive 
case only for finite temperatures, which are much higher 
than T cn and smaller than T c . 

It is also worthwhile to mention that the normal Fermi- 
liquid with the hard-core repulsion between fermions be- 
comes a p-wave [22[ or a d-wave [HI superfiuid at very 
low temperatures (on the mK scale) due to the Friedel os- 
cillations [24J of the particle-particle interaction potential 
caused by screening. This unconventional pairing should 
not affect our results for the proximity effect with the 
conventional s-wave superconductor as long as there is 
no symmetry-breaking mechanism. On the other hand 
if such a mechanism is involved (e.g. the spin-orbit cou- 
pling), the inhomogcneous s-wave order-parameter can 
generate a secondary order-parameter of another symme- 
try (2f| H^, [2?| . Finally, the order-parameter in a dirty 
normal metal should also exhibit some power-law decay 



at low temperatures, but the power could be different 
from 1/x of the clean metal |28j . 

In summary we have solved a long-standing problem of 
the low-temperature proximity effect in a clean normal 
metal with repulsive interaction between electrons. We 
have found the power-law decay of the order parameter 
with the characteristic proximity length strongly reduced 
with respect to the superconductor coherence length by 
the repulsive interaction in the normal metal. 
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